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Abstract 

The analysis of stochastic loss networks has long been of interest in computer and 
communications networks and is becoming important in the areas of service and 
information systems. In traditional settings, computing the well known Erlang formula 
for blocking probability in these systems becomes intractable for larger resource 
capacities. Using compound point processes to capture stochastic variability in the 
request process, we generalize existing models in this framework and derive simple 
asymptotic expressions for blocking probabilities. In addition, we extend our model to 
incorporate reserving resources in advance. Although asymptotic, our experiments show 
an excellent match between derived formulas and simulation results even for relatively 
small resource capacities and relatively large values of blocking probabilities. 
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1. Introduction 



The problem of satisfying a stream of customer (user) requirements from resources of 
finite capacities for some random processing time has long been present in many areas such 
as telephone and communication networks, inventory control (rental industry) and, recently, 
workforce management. For all of these applications, system dynamics can be described 
as follows. Requests for resources arrive according to some point process in time. If there 
are enough available (non-engaged) resources to satisfy their requirements at the moment of 
arrival, required resources are committed for some random time that represents their processing 
duration (holding time) after which they are released and become available to accommodate 
future requests. In the case of insufficient amount of available resources at the moment of its 
arrival, a request is lost. The previously described system is usually referred to as a loss net- 
work, and one of the commonly analyzed performance metrics is the blocking probability, i.e., 
probability that an incoming request is lost due to insufficient amount of available resources to 
satisfy its requirements. 

Loss networks with fixed resource requirements have been intensively analyzed in the con- 
text of circuit-switched networks. Let requests require resources of K < 00 different types for 
some random generally distributed processing time with finite mean. Furthermore, assume that 
requests belong to M different classes characterized by their resource requirements, processing 
durations, arrival rates. Then, assuming that requests of different types arrive according to 
mutually independent Poisson processes, by PASTA property ([24]), blocking probability £>; 
of an incoming request of type 1 < I < M is equal to the sum of probabilities of blocking 
states for I type request and is computed using the generalized Erlang formula (e.g., see [15]), 
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i.e., 



Bi = l-G(C)- 1 G(C-Ae ; ) 



where 



/ M 

g(c)= e n 



neS(C) 1=1 



Pi 



) 



and 



S(C) := {n G : An < C}, 



(1) 



where n = (m, . . . , n M ) and C = (Ci, . . . , Cr-). In the previous expressions Cfc, 1 < fc < 
if, is capacity of resource type k, A = [Ajy] is a if x M matrix, where Am represents the 
amount of resources of type 1 < k < K required by a request of type 1 < I < M, and pi, 
1 < I < M, represent traffic intensities of I type requests (computed as pi = Xi/pi, where 
A; is the arrival rate of I type requests and is the corresponding mean processing time). 
Furthermore, e; is a M dimensional vector with the Ith component equal to one and the rest 
equal to zero. In the case of a single resource type and a single request class with exponentially 
distributed processing times, blocking probability was first expressed by Erlang in 1917 (see 
[7]). Later on, it was shown that the Erlang formula holds under more general assumptions on 
call holding time distributions (see [20]) and in the case of Poisson arrivals with retrials (see 
[4]). It is noteworthy to point out the difference between the Erlang loss network and a queue 
with finite buffer. The two systems follow very different dynamics resulting in a different 
behavior and, therefore, their analysis (e.g., see [12] and [2]). 

It is easy to see that the cardinality of the state space 5(C) in (1) increases exponentially in 
the norm of vector C, i.e., |C| = Y^f=i | C* | • It is shown in [18] that the calculation of G(C) 
is a |JP-complete problem, which belongs to a class of problems that are at least as hard as 
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./VP-complete problems. To this end, many approximation techniques for evaluating blocking 
probabilities in large loss networks have been proposed. One of the most popular ones is 
known as Erlang fixed point method. The main idea of this approximation is to assume that 
deficiencies of different resource types happen independently. The application of the Erlang 
fixed point method can be traced back as early as 50's (e.g., see [23]). In [14], Kelly studied 
the performance of the Erlang fixed point method and established its relation to a nonlinear 
optimization problem. He also proved uniqueness of the fixed point and its asymptotic exact- 
ness when resource capacities and arrival rates grow with the same rate (see [15]). Some of the 
related practical aspects of Kelly's analysis were investigated in [22]. The Erlang fixed point 
method is further refined in [25]. There are also many other types of approximations such as 
recursive algorithm in [13], or unified approach based on large deviations for all (light, critical 
and heavy) traffic regimes in [8]. Overall, except from the bounds in [8], these methods make 
use of the structural properties of the Erlang formula and, hence, largely rely on the Poisson 
assumption for call arrivals. Another restriction of the above models is that the amount of 
resource requirements are assumed to be fixed; in fact, it is assumed that they are (0, 1) 
parameters in most of the cases considered. Meanwhile, we see in many applications that 
resource requirements could be highly variable and their distributions possibly long-tailed; for 
specific examples, see [10], [11] and [16]. Furthermore, more recently, loss networks models 
have been applied in the context of workforce management applications (see [19]), where 
requests behavior is even more volatile and extreme. 

In this paper, we analyze loss networks that have renewal arrivals and random resource 
requirements. In particular, we assume that request arrivals follow a compound renewal pro- 
cess, with the corresponding holding times being arbitrarily distributed with finite mean, in- 
dependent of each other and arrival points. In order to cope with variability in resource 
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requirements, we model them as subexponential random variables. We obtain a simple and 
explicit asymptotic expressions for blocking probabilities when capacities of resources grow. 
For the case of a single resource loss network, we show that the stationary blocking probability 
is approximately equal to the tail of the resource requirement distribution. In addition, we 
extend our results to allow advance reservations of resources. Finally, we investigate gen- 
eral (multiple resources and arbitrary topology) loss networks and show that the asymptotic 
blocking probability behaves as the tail of the heaviest-tailed resource requirement. Although 
asymptotic, our numerical experiments show an excellent accuracy of the derived formulas 
even for relatively small capacities and relatively large values of blocking probabilities, sug- 
gesting wide applicability of the obtained results. 

Our paper is organized as follows. In Section 2, we introduce our model in the context of a 
single resource type. Then, in Subsection 2.1, we state and prove our main result in Theorem 
1, while in Subsection 2.2, we extend it to the case of advance reservations. Further extension 
to the analysis of the stationary blocking probability in the case of general loss networks is 
stated and proved in Theorem 2 of Section 3. Our simulation experiments for some specific 
cases of arrival processes and resource requirements are presented in Section 4. Finally, we 
conclude our paper in Section 5. A discussion and the proof of existence of the stationary 
blocking probability is presented in the Appendix. 

2. Systems with one resource type 

Let requests for resources from a common resource pool of capacity C < oo arrive at time 
points {t„, — oo < n < oo} that represent a renewal process with rate < A < oo, i.e., 
E[r„ — r n _i] = 1/A. At each point r„, B n amount of resources is requested. If available 
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capacity is less than B n , this request is rejected (blocked); otherwise, it is accepted and B n 
amount of resources will be occupied for the length of time 6 n . Sequences {B n } and {6 n } of 
i.i.d. random variables (r.v.) are assumed to be mutually independent and independent of the 
arrival points {t„}; furthermore E6 n < oo for all n. Let B and 9 denote random variables that 
represent {B n }, {8 n }, i.e., F[B > x] = ¥[B n > x],F[0 > y] = F[6 n > y], for any n G Z, 
x > and y > 0. 

In this paper, we assume that B is a subexponential random variable, defined as follows 
(e.g., see [9]): 

Definition 1. Let {Xi} be a sequence of positive i.i.d. random variables with distribution 
function F such that F(x) < 1 for all x > 0. Denote by F(x) = 1 - F(x), x > 1, the tail 
of F and by F n * = 1 — F n * (x) = Ppfi + • • • + X n > x] the tail of the n-fold convolution 
of F. F is subexponential distribution function, denoted as F e S, if one of the following 
equivalent conditions holds: 

• Hindoo F F ff = n for some (all) n > 2, 



lim T 



oo m^ Xl tV" +X v i = 1 for some (all) n > 2. 

^° IP[max(Xi ,...,X n )>xJ v y — 



For a brief introduction to subexponential distributions the reader is referred to a recent 
survey [9]. This class of distributions is fairly large and well known examples include regularly 
varying (in particular Pareto), some Weibull, log-normal and "almost" exponential distribu- 
tions. 

Next, let Afn°^ be the set of indices i < n of resource requirements that arrive prior to r„, 
are accepted, and are still active by time t„. Furthermore, let = \Afi C ^ | be a cardinality 
of set Mn ■ Thus, the total amount of resources Q n that an arrival at time t„ finds engaged 
can be expressed as = J2 ie j^( c "> B%- 
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P[Qf) + B n > C] 



(2) 



for large C. It can be shown that for the model introduced above there exists a unique stationary 

(c) 

distribution for Q n and, therefore, the quantity in (2) is well defined. The proof of this 

( c) 

result is based on constructing a Markov chain with general state space, of which Q n is a 

functional. Then, by using a discrete version of Theorem 1 from [20], we show that there exists 

( c) 

a unique stationary distribution for the constructed Markov chain (and, therefore, Q„ ') which 
is ergodic. Since this proof is not the main focus of this paper, we present it in the Appendix. 

In this paper we use the following standard notation. For any two real functions a(t) and 
b(t) and fixed t G RU{oo}, let a(t) ~ b(t) ast—>t denote lim t ^ to [a(t)/b(t)] = 1. 

2.1. Blocking probability in a system with one resource type 

In this section we estimate the stationary blocking probability P[Ql C ' + B n > C] in a loss 
network with a single resource pool when its capacity C grows large. 

Theorem 1. Let {B n , —oo<n< oo} be a sequence of subexponential random variables 
with finite mean. Then, the stationary blocking probability satisfies 



Proof: First, observe that a request will be lost if it requires more than the total capacity C 
and, therefore, 



P[QP + B n > C] - V[B >C] as C -» oo. 



(3) 



P[Qi c) + B n >C}> F[B > C] for all C> 0. 



(4) 
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In order to prove the asymptotic upper bound for P[Ql C ' + B n > C], we start by condi- 
tioning on the size of B n as 



P[Ql C) +B n >C]= P[Ql c > + B n > C, B n > C] + P[Ql c ) + B n > C, B n < C] 



= h+h- 



(5) 



Note that h is upper bounded by F[B > C}. Next, we prove that I 2 = o(¥[B > C}) as 
C — > oo. In view of the definition of A/"i°' from above, 



(6) 



Observe that for i G A/"„ , BiS are mutually dependent which makes direct analysis of the 
expression in (6) complex. For that reason, we sample the original process of arrivals at points 
Tj at which the requested amount of resources 5, is smaller or equal to C and observe another 
system of unlimited capacity with the sampled arrivals. Let N n , s be a set of request indices 
i < n that belong to the sampled process and are still active at time r„, i.e., 

Af s ,n = {i < n\Bi <C,9 l > r n - Ti}. 

Note that the sampled process is renewal as well with rate XF[B < C]/P[B > C] and 
that resource requirements B i7 i e Af 8 ,n> are mutually independent. Furthermore, since 
Nn C A/ S>n , we can upper bound 7 2 in (6) by the probability that the total amount of 
required resources in a new system exceeds capacity C, i.e., 



h < 



Bi + B n > C, B n < C 



(7) 



Now, in view of the results derived in [6] for every integer n and i.i.d. subexponential random 
variables B u ..., B n , P[ElLi B i > C ) ~ P[max(Bi, B 2 , . . . , B n ) > C] as C -» oo, 
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implying asymptotic relation 

n 

P ^B t > C,B t <C for every 1 < i < 



o{P[B >C]) as C -> oo. 



In order to show that n can be replaced by 7V Si „ in the above inequality, we need to integrate it 
with respect to the density of N Sin , i.e., 



V" Bi > C,Bi < C for every i E W s ,„ U {n} 

'fc+1 



^P[JV s ,„ = fc]P 



k=0 



^2 Bi > C, Bi < C for every i = 1, . . . , k + 1 



i=i 



Note that on the left hand side of the previous equation index i can take negative values. Next, 
due to the lemma stated by Kesten (see Lemma 7, pp.149 of [3]), for any e > there exists a 
positive constant K (e) such that 

P[Ei=i Bi > C, B t < C for every 1 < i < k] P[£ -=i B t > C] 



>C] 



< 



>C] 



<K(e)(l + e) k , 



for any integer k and all capacity values C < oo. Then, since the probability generating 
function Kz Ns - n is finite for any z £ C (see Theorem 1 in [21] and Theorem 5 in [17] for 
the detailed proof), we have J2kLo ^Ws,n = + e) fc < oo. Therefore, by applying the 
dominated convergence theorem, we conclude that 

p [E ie w 3 „ Bi + B n > C, B, < C for every i e Af s , n U {n} 
lim — — 



= iim y 



P[B > C] 

P[N Stn = k}¥ \j2i=i Bi > C, Bi < C for every 1 < i < k + 1 



fe=0 



¥[B > C] 



0, 



(8) 



which in conjunction with (5) and (4), completes the proof of this theorem. 



O 
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Remark: It may appear surprising that the performance of the loss network from above does 
not depend on engagement durations, as long as they have finite mean. In addition, the result is 
quite general and provides the asymptotic result for a large (subexponential) class of possible 
resource requirement distributions. 

2.2. Advance reservations 

Using the result of Theorem 1 and observations from the previous remark, we extend the 
loss networks model to allow requests to become effective with some delay with respect to 
the moments of their arrivals. In particular, a request that arrives at time r„ and requires B n 
amount of resources for some random time 8 n starting from the moment r„ + D n is accepted if 
previously admitted resource requirements allow that; otherwise, it is rejected. In other words, 
a request arriving at t„ is lost if at any moment of time in interval (r„ + D n , r„ + D n + 9 n ) 
the total amount of active requirements requested prior to r„ exceeds C — B n . First, note 
that B n > C implies the loss of nth request and, therefore, it is straightforward to conclude 
that the blocking probability in the system with advance reservations can be lower bounded by 
F[B > C}. 

Next, we discuss the idea behind proving the upper bound on the blocking probabilities. 
By applying sample path arguments one can show that, at any moment of time, the amount of 
active resources in the previously described system with advance reservations can be bounded 
from above by the amount of active resources in another system of unlimited capacity, without 
advance reservations, with resource holding times D n + 9 n for every n, and with requests 
for resources being sampled from the original process {B n } whenever the corresponding 
requirements are less or equal to C. Equivalently, the blocking probability in the system with 
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^ Bi + B n > C 

ieAf^J (B+D) 



where N$(0 + D) is a set of request indices i < n that are active at time r n , whose 
requirements are less or equal to C and holding times last throughout the interval (r, , n + 
Di + Oi), assuming that there is an unlimited resource capacity. 

Finally, by using the previous discussion, the properties of {B n }, {9 n } and {t„} as intro- 
duced at the beginning of this section, assuming that reservation times {D n }, ED n < oo, are 
i.i.d. and independent from {B n }, {9 n } and {r„}, and applying the identical arguments as in 
the proof of Theorem 1, we obtain the following result: 

Corollary 1. The blocking probability in the system with advance reservations approaches 
F[B > C] as C -> oo. 



3. Acquiring resources of different types (loss networks case) 

Assume that there are K e N resource types with capacities C\, . . . , Ck- Again, requests 
arrive at {r„, — oo < n < oo}, which represent a renewal process with rate < A = 1/E[ri — 
r ] < oo. There are M < oo request types and, given an arrival, the request is of type I, 
1 < I < M, with probability pi, pi + ■ ■ ■ + pm = 1. independent from {t„}. We will 
use random variables J„ 6 {1,2,..., M} to denote the type of the request arriving at t„. 
Furthermore, let B^ Jn ' 1 \ . . . , Bn 7 "^ represent amounts of required resources of each type 
at time r„ and let 9n"\ E0ri J "^ < oo, be the corresponding random duration. We assume 
that sequences {(Bn n,1 \ . . . , £?i J "' K ')}, {9n Jn ^} are mutually independent and independent 
from {t„}. Given the event { J n — I}, resource requirements Bn' l \ 1 < i < K, are mutually 
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independent nonnegative random variables drawn from distributions Fij, 1 < i < K; if a 

(l i) 

request does not require resources of type i then B. n = a.s., — oo < n < oo. Only if 
there is enough capacity available, the request arriving at time t„ will be accepted and all of 
the engaged resources will be occupied for the duration of #i J "'; otherwise, the request is 
rejected. 

Our goal is to estimate the blocking probability in a system described above. Define 
Qn \ ■ • • , to be amounts of resources of each type that a request arriving at time t„ 

(i) 

finds engaged. Note that Q„ , 1 < i < K, are mutually dependent and, as pointed out in the 
Introduction, it is hard to compute the blocking probability of this system explicitly. Using 
analogous arguments as in the case of a single resource type (see the Appendix), one can show 

(i) 

that the stationary distribution of Q\ ! , 1 < i < K, exists. Probability that the request arriving 
at time r„ is blocked equals to 

P[Ui< i < Jf {g«+S( J »' < ) >d}], (9) 

and our goal again is to estimate its value as minj d grows large. 

Asymptotic estimates derived in this section hold under the following assumption: 

Assumptions: For each resource type 1 < i < K, let d and Hi be two disjoint sets of request 
types (|£j \JHi\ = M) satisfying: 

• Assume that there exists at least one resource type that is accessed by subexponentially 
distributed resource requirements, which implies \Hi\ > for some 1 < i < K; 

• For every I e Hi ^ 0, there exists a subexponential distribution Fj e S such that 

Fi,i(x) ~ ci yi F l (x) as x -> oo with c iti > 0; 
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• There exists a subexponential random variable L e S that satisfies 

P\L > x}> max PLB^"^ > x\ J n = 1} for all x > 0, 

l<i<K,leCi 

and F[L > x] = o{F,{x)) as x -> oo for all i e {j\Hj ^ 0}. 

Remark: In the preceding assumptions, we require the resource requirement distributions to 
be asymptotically comparable. For each 1 < i < K,Hi contains tail dominant subexponential 
distributions that are asymptotically proportional to each other. On the other hand, the only 
assumption imposed on the distributions in d, 1 < i < K, is that there is a subexponential 
tail that asymptotically dominates them. This asymptotic tail comparability is necessary for 
our main result to hold. In particular, these conditions are extensively used in (16) - (21) of the 
proof of Theorem 2. 

Next, we prove the following lemma that investigates summations of random variables with 
different tail distributions. 

Lemma 1. Suppose that Xi, 1 < i < n, are independent random variables with correspond- 
ing tail distributions Fi(x), 1 < i < n. If there exists F G S such that Fi(x) ~ CiF(x) as 
x — > oo with Ci > 0, 1 < i < n, and X^Li Ci ^ ^» t ^ len the following asymptotic relation 
holds: 

n 

= o(F(x)) as x -> oo. (10) 



Xi > x, Xi < x, 1 < i < n 

.i=l 



Proof: Note that 



n 




^Xi> x 


= P 







Xj > x, Xi < x, 1 < i < 



+ . 



J2x i>X ,UU{Xi>x} 



i=i 



Then, the previous expression, U" =1 {Xi > x} C > x }> independence of JQs, as 

well as Lemmas 4.2 and 4.5 of [1], imply (10). O 
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First, we estimate the asymptotic lower bound for the expression in (9). By using our model 
assumptions, {Bn n ^ > d} C {Qn^ + Bn"^ > C} and independence, we obtain 

K 

P[Ui<i<*{Q« + B^ > d}] > P[Ui< i < Jf {B( J »' < ) > Q}] ~ E PiFiACi), 

»=i lew, 

(11) 

as mini C — ► oo. 

Next, we estimate the asymptotic upper bound for the expression in (9). Using the union 
bound yields 



K 

n^<i<K{Q { : ] + 4 Jn,<) > ct}] < Y, ¥ iQn + B n n ' l) > d]- ( 12 ) 

i=l 



Similarly as in (7) of Theorem 1, for each resource 1 < i < K, 



nQ^+B^ >d] < 



E E B f i] + E E Bf) + ^-.o>c i 



(13) 



where Af^rf'K 1 < i < M, are sets of indices j < n defined as 

i} = {.7 < n| J, = Z, i} < Cf, (if > r n Tj }. 

In the previous expressions we bounded the amount of allocated resources that are active 
at time r„ by the corresponding quantity in another system of infinite capacity where the 
corresponding request process is sampled from the original {£?i' 7 ™'^}, 1 < i < K, whenever 
the corresponding requirements are less than or equal to C 1 < i < K. 

In the rest of the proof, we derive an asymptotic estimate for the expression in (13). After 
conditioning on {N& Ct) = m, . . . , N { s %' Ci) = n M } (N { s l f l) = \^% C ' ] \, 1 < I < M), we 
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obtain 



P[Q^+Bk JnA >Ci] 

< J2 nNi]n Ci) =nu...,N^)=n M ] 

0<ni,...,riM<oo 



E E B ii) + E E B uf + B ^' l) > C " B U) <C*,l<3<m,l<l<M 
le-Hi j=i led j=i 

(14) 



where = B< M) , k S Afsln'*\ j — 1, . . . , ni, are independent replicas of requests in 

Ns]n^- Next, after conditioning on { J n = m}, m = l,...,M, and then on B^' 1 ^ being 
smaller or larger than d, we can further upper bound the conditional blocking probability in 
(14) as 



E E B W + E E B W + B n n ' l) > B W <Ci,l<j<n h l<l<M 
ien t j=i led 3=1 



< 



M 

E^ 

m— 1 



E E B ti) + E E B u) ] + B n m ' l) > c *> B U ] ^ Ci - 1 ^ j — n i< i — i — m, 4 m,<) < Ci 

ieHij=i led 3=1 



M 



(15) 



Thus, the probabilities in the first term on the right hand side of the previous expression can be 
expressed in the form 



E E B U? + E E B U) ] > C " B U? <Cul<j<n' h l<l<M 
leHi j=i led 3=i 



, (16) 



where n[ = ni for I ^ m and n\ = n\ + 1 for I = m. 

Next, in order to estimate the asymptotic upper bound of the term in (16), Assumptions 
enable us to distinguish between two cases: (i) Hi — or J2ieH - n 'i ~ 0> an( ^ ^* ^ ^ an d 
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(i): If ?U = or X/ie'Hi n l ~~ ^' we h ave th a t probability in (15) can be upper bounded as 



< 



where in the inequality above we used Assumptions and introduced L^hf 1 to be independent 

(l i) 

r.v.s equal in distribution to L. Hence, since ; are subexponential, we obtain 



lim 



<E< 



(17) 



(ii): If Wj and X^ eW . n 'i > ®> usm 8 Assumptions and Lemma 1, we derive the following 
asymptotic upper bound 



EE 5 W + EE S 8> ) >Cl,Bgf <Cul<j<n' u l<l<M 

1 /- u a — i i /- r . A — i 



lent j=i 



led o=i 



o(Fi(Ci)), 
(18) 



as d — ► oo. 

Thus, in (16)-(18) we obtained upper bounds and their asymptotic estimates for the con- 
ditional blocking probabilities in the first term of (15) that hold for any finite nonnegative 
integers m, . . . , um- Thus, in view of (14), in order to estimate an asymptotic upper bound of 
P[<2n^ + Bn Jn ' 1 ^ > d], we need to integrate probabilities in (16) with respect to densities 

(I C ) 

of r.v.s Ng y k , I = 1,...,M. In this regard, note that in the case where Hi ^= 0, by 
Assumptions, the term in (16) can be upper bounded as 



E E B Sf + E E B !Jf > B (i) } < a, i < j < «{. i < / < m 

leHi j=i led j=i 



< 



n l n l 

EE B 8) <) + EE L S) )>Ci 



(19) 



ien z j=i led 3=1 

(l i) 

where, as before, are independent r.v.s equal in distribution to L. Furthermore, since 

F[L > x] — o(Fi(x)) as x — > oo, there exists a large enough finite integer if such that 



Asymptotic Blocking Probabilities in the Loss Networks with Subexponential Demands 



17 



F[L > x] < HFi(x) for all x > 0. Therefore, for any x > 0, one can write 



F[L >x}< HFi{x) = P Ui< r <H{S« > x} 



< 



H 



E^ 



> X 



(20) 



where Br, 1 < r < H, are independent r. v. s having cumulative distribution function Fj. 
Now, in view of (20), each of random variables L y J^' in (19) can be stochastically upper 
bounded by a random variable that is equal in distribution to J2r=i • Thus, if we introduce 
Yj, j > 1> to be independent r.v.s equal in distribution to ^ r =i ^ > we obtain 



EE B S) i) + EE L S) )>Ci 

ien z j=i led j=i 



< 



EE*Sf + E *s><* 



which in conjunction with point (b) of Lemma 4.2 in [1] implies that for any e > there exist 



a finite constant K e such that 



EE B S) i) + EE £ S) )>c * 

ieHij=i ieCij=i 



< 



EE*Sf + E 

ie-Hi j=i i=i 



< # e (l + e) Ele "i " i+E '«. n ' l F t (Ci)i (21) 

for any Cj < oo. Similarly, in cases where Wj = 0, we could apply the stochastic dominance 
^ ^(j)^' ' e wnere ^(j)*'' are ' as before, independent subexponential random 
variables equal in distribution to L. Then, by Kesten's lemma (see Lemma 7 on page 149 
of [3]), the analogous bound to the one in (21) follows. 

Finally, since (21) bounds uniformly probabilities in (16) for all d < oo and n\, 1 < I < 
M, in conjunction with (15), (14), N^ i] < iV^' oo) a.s. and existence of Ez N ^°° } for all 
z G C, 1 < I < M, (see Theorem 1 in [21] and Theorem 5 in [17]), one can apply the 
dominated convergence theorem and conclude 

lim ^n^ > Q] 

Next, by adding asymptotic estimates for all 1 < i < K, in conjunction with (1 1), we complete 
the proof of the following result: 
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Theorem 2. For the request model introduced in this section, under the conditions imposed 
by Assumptions, the stationary blocking probability for general loss networks satisfies 

K 

lP[Ui<i<_ff {Q« + > d}} PlClAid) as minC, - oo. 

1=1 l£Hi 

4. Numerical examples 

In this section, with two simulation experiments, we demonstrate the accuracy of our 
asymptotic formulas, proved in Theorems 1 and 2. Our goal is to show that even though 
our results are asymptotic, the derived estimates match experiments with high accuracy even 
for systems with finite support demand distributions and moderately large capacities. 

In each experiment, in order for the system to reach stationarity, we let the first 10 8 arrivals 
to be a warm-up time. By repeating many experiments, we observe that longer warm-up times 
do not lead to improved results. Then, we count the number of blocked requests among next 
10 9 arrivals. In both of the experiments below, measurements are conducted for capacities 
C = 500 + lOOj, < 3 < 9, where the starting value of C = 500 is set to be slightly larger 
than the effective systems load AE[0 n ]E[B„]. Simulation results are presented by symbol "o" 
in Figures 1 and 2, while our approximations, estimates obtained in Theorems 1 and 2, are the 
solid lines on the same figures. Note that in order to emphasize the difference and to observe a 
range of blocking probabilities we are trying to estimate, we present base 10 logarithm of the 
obtained values. 

Example 1 Consider the case of a single resource type of capacity C. Let requests for resources 
arrive at Poisson time points with rate A = 1. In addition, we assume that engagement 
durations are exponentially distributed with mean 1/fi = 1. Next, let request requirements 
B n be drawn from a finite support distribution, where F[B n = i] = 1 < i < 1999, and 
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¥[B n = 2000] = 1-P[B„ < 2000] (power law distribution). Effective load in this example is 
AE[0 n ]E[_B n ] » 485.8. Experimental results are presented in Figure 1. Even though we start 
measuring rejections at capacities that are slightly larger than the mean requirement value, our 
approximation P[B„ > C] is very close to experimental results. In particular, the relative 
approximation error is less than 1% for C = 500, and for capacity values larger or equal to 
C = 1400 this error is less than 0.3%. 

-0.605 

-0.61 

g -0.615 
h 

a 
o 

o 

2 -0.625 
-0.63 
-0.635 

500 600 700 800 900 1000 1100 1200 1300 1400 
resource capacity C 

FIGURE 1: Illustration for Example 1 

Example 2 In this example, we consider the case of two resource and two request types. 
Furthermore, we assume that resource capacities are the same C = C% = C%. The frequencies 
of requests of types 1 and 2 are p\ = 0.3 and p2 — 0.7 respectively. Assume that the arrival 
points are separated by a fixed, unit length of time, i.e., r n — r n _i = 1 for all n. Type 1 request 
durations satisfy 9^ ~ exp{4) and type 2 request holding times are drawn from the uniform 

(2) 

distribution on [0,40], i.e., 0\ ~ Unif([0, 40]). Resource requirements corresponding to 
engagements of type 1 are distributed as P^ 1 - 1 = 1] = 0.8, P[.B M = i] = 0.15e"A 
2 < i < 1999 and PfB 1 ' 1 = 2000] = 1 - P^ 1 ' 1 < 2000] for the type 1 resources, and 
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F[B 1 ' 2 = 50] = 1 for type 2 resources. Requests of type 2 require resources according to 
¥[B 2 ^ =i]= georri^il - geom), 1 < i < 1999, V[B 2 ^ = 2000] = 1 - P^ 2 - 1 < 2000], 
where geom = 0.6 for resources of type 1, and F[B 2 ' 2 = i] = £4, 1 < i < 1999, P[5 2 < 2 = 
2000] = 1 — P[£? 2 ' 2 < 2000] for type 2 resources. Our asymptotic results suggest that the 
blocking probability should be characterized by the heaviest tailed demand distributions. The 
results of this experiment are presented in Figure 2. As in the previous case, we obtain a very 
accurate agreement between our approximation and the simulation. The relative approximation 
error in this case does not exceed 2% and is getting smaller as resource capacities grow. 

-0.71 
-0.72 
-0.73 
1 -0.74 



-0.77 
-0.78 
-0.79 

500 600 700 800 900 1000 1100 1200 1300 1400 
resource capacity C 

FIGURE 2: Illustration for Example 2 

Remark: (i) We would like to point out that the accuracy of experimental results directly 
depends on the approximation errors (7) and (18), depending on the simulated scenarios. These 
errors highly depend on the tail properties of the resource requirements distributions. More 
specifically, under fairly general assumptions, the heavier the dominant tail of the resource 
requirement distribution is, the smaller would be the relative approximation error. For detailed 
explanations, a reader is referred to Section 1.3.2 of [5]. (ii) Note that our main results 



asymptotic formula 

o simulation 



A 

o 

o 



o 
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estimate the stationary blocking probability and, as we commented earlier, are indifferent 
to distributional properties of holding times. For that reason, as long as one can claim that 
the measurements are conducted in stationarity, the transience should not affect experimental 
results. 



5. Concluding remarks 

In this paper, we consider loss networks with reusable resources and finite resource ca- 
pacities and estimate the probability that a request is rejected due to insufficient amount of 
resources at points of their arrivals. Assuming a renewal process of request arrivals, subex- 
ponential resource requirements and generally distributed activity durations, we show that the 
asymptotic blocking probability for a wide class of analyzed systems can be fully estimated 
using resource requirement distribution, independent from other system's properties. In par- 
ticular, we show that the blocking probability behaves as the asymptotically dominant tail of 
the resource requirement distribution. 

The model we study can be applied to a wide range of applications. Historically, loss 
networks (in particular, Erlang loss networks) are widely used for modeling communication 
networks. Later, through the development of new services applications such as workforce man- 
agement with similar modeling properties, the importance of accurately estimating blocking 
probabilities of general loss networks has become significant. In this regard, we investigate 
loss networks with various request types and possibly highly variable random amounts of 
required resources. In addition, we research the possibility of incorporating random advance 
reservations for incoming requests. These results should be of great interest to an emerging 
research community. Although our results are intended mainly for qualitative purposes, nu- 
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merical examples demonstrate an excellent match between derived formulas and simulated 
systems performance, hence strongly suggesting their application. 



Appendix 

In this section we prove the existence of the stationary blocking probabilities in (2). Using 
the model description from Section 2, we observe the system at the moments of request 
arrivals. Then, we define a discrete time process X n = (Ni C \Bi, Ei, i = 1, 2, • • • , Nq C J), 

where EiS represent times that elapsed in processing requests in the system by time r„; 

(c) 

furthermore, Nq J is the number of active requests at the moment of nth arrival. Note that 
{X n }„> is a discrete time Markov chain with state space £1 = No x ioc x l^, where 
denotes the Banach space of the infinite sequence of real numbers equipped with the supreme 
norm; let uu e £1 denote the state with no active requests. We start observing the system at the 
moment t of Oth arrival and denote the initial state by X Q = (N°,B°,E°, i = 1,...,N°) 
drawn from some arbitrary distribution P , where E£?° < oo, E6*° < oo. Next, define T to 
be the Borel field of Q, and let P n (xo, A), x E ft, A e T, represent a transition probability 
of the Markov chain X n into set A in time n, starting from state xq. Let P n be the probability 
distribution of X n . 

Now, in order to prove the existence of a unique stationary distribution for the Markov 
chain {X„}, we use a discrete version of Theorem 1 in [20], which we state next for reasons 
of completeness. 

Theorem: A Markov chain homogeneous in time has a unique stationary distribution which 
is ergodic if, for any e > 0, there exists a measurable set S, a probability distribution R in Jl, 
and ni >, k > 0, K > such that 
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• kR{A) < P ni (x, A) for all points x E S and measurable sets A C S; for any initial 
distribution P there exists n such that for any n > n 

• P n (S) > 1-e, 

• P n {A) < KR(A) + e for all measurable sets Ac S. 

Proof: The proof follows identical arguments as in [20] translated into discrete setting. 

Next, we need to show that Theorem 1 holds for the process investigated in this paper; 
in particular, we will consider a common resource pool case. The proof follows the similar 
reasoning as in Theorems 4 and 5 of [20]. 

Define set S(ip, (3, 8) as 

Sty, /?, 8) 4 {N$ < V, < Bi < (3, < E t < 8, i e Af$} 

for some positive finite constants tp, (3, 8. 

Now, we show that for any e > 0, there exists S(ip, (3,6) C ft such that for any initial 
distribution P there exists n such that for all n > n 

P n (S(iP,/3,6))>l-e. (22) 

Note that 

P n {S{4>,f3,8)) < F[U ieA ,(c,{^ > 8},N<g < +P[U. £A/ ,o){ J B i > /3},N<g < 4>)+P[N^ > $\ 

< >s}+ >p}+ mi c j + <„ > ^\ , w 

(C) 

where Na,n represents the number of active requests at r„ that originated from n arrivals at 
r , . . . t„_i, and the rest of active requests at r n , N® n — Nq C J — N^n are those that were 
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active at the initial point r and are still processed at the moment of nth arrival. Next, since 



n 2 



N (C) > t 



JV 0,n > 2 



< 



< 



N" 



> Tn - To] > 



i=l 



p [ N ° = m \ v 



m=0 



2l[(9?>(l-c 1 )nE[ri-7 l ,]]> 



+ P[r„-r < (1 - ei)nE[n - T ]], 



(24) 



where in the previous inequalities < e\ < 1 is an arbitrary constant and we used Nn°°^ > 
Na°n a.s., where is defined as in the proof of Theorem 1. 

Now, we prove that there exists ip — ipo large enough such that (24) is bounded by e/3. By 
definition of Ni 00 ^ in Section 2 and Little's formula, EiV,! 00 ^ < oo and, therefore, 



lim P 



AT(°°) > t 



0, 



(25) 



uniformly for all n > 0. Next, note that 1 [0? > (1 - ei)nE[n - r ]] < 1[0° > (l-ei)E[n- 
tq]] a.s., and that for any fixed m, 



X;i[fl?>(l-ei)nE[ri-7D]]> 



i=l 



< 



^l[^ >(l-e 1 )E[r 1 -ro]]> 



| as V — * oo, 



which by the monotone convergence theorem implies that the second term in (24) satisfies 



lim V V[N a = m]F 

m=0 



£l[(??>(l- ei )nE[r 1 -7 l) ]]>| 



0, (26) 



uniformly for all n > 0. Finally, by the Weak Law of Large Numbers, for all n large enough, 



\tu - to < (1 - ei)nE[n - r ]] < e/9. 



(27) 



Thus, the previous conclusion in conjunction with (26), (25) and (24) implies that for an 
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arbitrary < e < 1, there exist n < oo and ip < oo large enough such that for all n > n 



< | and P t Ar ° = m ] p 

m=0 



J2l[0°>(l-e 1 )nE[T 1 -To}}> 



<1. 

- 9 



(28) 

Now, since EE-i < oo, E6*j < oo, there exist P , S , such that 

F[B, > A)] < and P[0 4 > *„] < 

Thus, the previous expressions in conjunction with (28), (27) and (23) imply that for all large 
enough n > n inequality (22) holds for a chosen set S(ipo, Po, ^o)- 

Next, we show that there exists m > and k > such that for all points x 6 S(i(>o, Po, So) 
and measurable sets A C S(ipo, Po, So), the following inequality holds 

P ni (x,A)>kR(A). (29) 

Let Fg(u) denote a cumulative distribution function of a random duration 6, i.e., F[6 < u]. 
Furthermore, select a small positive number r\ such that for some chosen A > 5 , Fg(A) — 
Fe(5o) — i] > 0. Next, for any n\ 

P ni (x,A)>P 1 (x,Lu )P n2 (uj ,A), (30) 

where n 2 = n x - 1. Let x = (m,bi, . . . ,b m , a, . . . ,e m ) € S(ip , Pa, S Q ). Then, 

P\{x, ivo) > — t~o > A, all m requests depart in (to, n)] 

J a fJi 1-F e (ei) 
where F a (u) represents cumulative inter arrival distribution of a renewal process {r„}, i.e., 
F a (u) = P[n — r <u]. Now, by applying lower bound 

F e {u + e i )-F e {e i ) 



l-F e ( ei ) 



> F e (A) - F e (6o) = v 



26 Lu, Radovanovic 

in (31) we obtain 

Pi(x,wo) > v m P[ri - ro > A] > ^(1 - F a (A)). (32) 
Next, we derive a lower bound for P n2 (uj , A) for some n 2 large enough such that 

P[t„ 2 -t ><5]>1-|. (33) 

Note that the condition imposed on n 2 in (33) is possible due to the Weak Law of Large 
Numbers, since for any e > and all n 2 large enough with Sq < (1 — e)E[r„ 2 — t ], 

P[r„ 2 - to > So] > P[r„ 2 - r > (1 - e)E[r„ 2 - r ]] > 1 - |, 

Next, pick any x' = (m', e' 1; . . . , e^,, 6' 1; . . . , b' m ,) G A where, without loss of generality, we 
assume that e[ > e' 2 > ■ ■ ■ > e' m ,. Define x' + dx' = (to', e[ + de[, . . . , e' m , + de^,,^ + 
db[, . . . , b' m , + db' m ,) where de[, • • • , de! m ,,db' x , ■ ■ ■ , db' m , are infinitesimal elements. Then, 
the transition probability into state (x',x' + dx') starting from uj can be bounded by the 
probability of the event that there are exactly to' arrivals prior to r„ 2 whose arrival times are 
determined by e' l7 • • • , e' m ,, whose resource requirements are in (b' 1 ,b' 1 +db' 1 ), ■ ■ ■ , (b' m ,,b' m ,+ 
db' m ,), and where the rest of n 2 — to' arrivals are rejected since their requirements exceed 
capacity C. Therefore, 

P n2 {wo,{x\x' + dx')) > |ll P ^ > e 'M B i e M + ^)]J p [ B i >Cr- m ' 

{'i2 — 1 ri2 — 1 

]T Yj G (ei - e' 2) ei - e' 2 + dei), ■ • • ^ Y j G (e' ro „ e^, + de' m ,) 
3=h 3=i m i 



«2 — 1 n 2 - 1 



U ^ E ^ G ( e 'l ~ e 2' 6 'l ~ 6 2 + de 'l)' ■ • ■ ' ^ G ( e m'' e m' + de '™>) 
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where / = {0 < i\ < ii < ■ ■ ■ < i m i < n 2 — 1} and Yj are i.i.d. random variables equal 
in distribution to inter-arrival times of the renewal process {r„}, i.e., Yj = Tj + \ — tj. Now 
denote 

Jlnejye'^iBjGib^b'j+db'j)} 

i'2 — \ n.2 — 1 



X P 



(34) 



and define probability distribution 



R(A)±v[ r(m , ,ei,...,e , m „6i,...,C). (35) 



where V is a normalization constant. Note that R(A) is well-defined since 

r(m',ei, . . . ,4,&i, . . .6^,) 

>o 



s / / 

< P[iv£ } < +P[^l oo) > Vo] < oo. 

The previous inequalities, in conjunction with (34), (32) and (30) imply that 

P n2+1 (x,A) > r]^°(l - F a (A))P[Bi > C^V^RiA). (36) 

Finally, it is left to show that there exists K > such that for every initial distribution Pq, 
for all n large and for any measurable set A C S{ipQ, /3q, Sq) 



P n (A) <KR(A) + e. 
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By (33), for all n > n 2 



P n {A) < P[X n e A, r n - T > 5q] + P[r„ - t < S ] 
< F[X n g A, t„ - T > <5 ] + e/2 



< 



x'eA 



J =1 



I 



12-1 



n 2 — l 



r(m / ,ei,...,ej n ,,6i,...,6^,) + e, 



where the second inequality follows from the fact that requests that are active at t„ must occur 

in the previous <5 length of time that are captured in n 2 renewal intervals [t„_„ 2 , t„_„ 2+ i), . . . , [t„_i,t„) 

with significant probability (greater than 1 — e/2). Thus, after applying definition (35), we 

obtain that for all n large 



P n (A)<V- 1 R(A) + e, 



which, in conjunction with (36) and (22), implies that the process X n satisfies conditions of 

the theorem stated at the beginning of this section. Thus, there exists a unique stationary dis- 

(c) 

tribution for the Markov chain X n . Therefore, since Q n defined in Section 2 is a functional 
of the process X n , it has a unique stationary distribution as well implying the existence of the 
stationary blocking probability. 

O 



Asymptotic Blocking Probabilities in the Loss Networks with Subexponential Demands 



29 



Acknowledgments 



The authors would like to thank Prof. Predrag Jelenkovic for valuable suggestions related 
to the possible generalizations of this work. 

References 

[1] ASMUSSEN, S., HENRIKSEN, L. F. AND KLUPPELBERG, C. (1994). Large claims approximations 
for risk processes in a Markovian environment. Stochastic Processes and their Applications 54, 29- 
43. 

[2] ASMUSSEN, S. AND PlHLGARD, M. Loss rates for Levy processes with two reflecting barriers. 
Mathematics of OR. to appear. 

[3] ATHREYA, K. B. AND Ney, P. E. (1972). Branching Processes. Springer- Verlag. 

[4] BONALD, T. (2006). The Erlang model with non-Poisson call arrivals. Procedings of 
ACM/Sigmetric and Performance Conference, 276-286. 

[5] EMBRECHTS, P., K. C. AND MlKOSCH, T. (1997). Modelling Extremal Events for Insurance and 
Finance. Springer. 

[6] EMBRECHTS, P. AND GOLDIE, C. M. (1980). On closure and factorization properties of 
subexponential and related distributions. J. Austral. Math. Soc. Series(A) 29, 243-256. 

[7] ERLANG, A. K. (1917). Solution of some problems in the theory of probabilities of significance in 
automatic telephone exchanges. Elektrotkeknikeren 13, 5-13. 

[8] GAZDZICKI, P., LAMBADARIS, I. AND MAZUMDAR, R. R. (1993). Blocking probabilities for 
large multirate Erlang loss systems. Advances in Applied Probability 25, 997-1009. 

[9] GOLDIE, C. M. AND KLUPPELBERG, C. (1998). Subexponential distributions. In A Practical 
Guide to Heavy Tails: Statistical Techniques for Analysing Heavy Tailed Distributions, ed. M. T. 
R. Adler, R. Feldman. Birkhauser, Boston pp. 435^-59. 

[10] HEYMAN, D. P. AND LAKSHMAN, T. V. (1996). Source models for VBR broadcast-video traffic. 
IEEE/ACM Transactions on Networking 4, 40-48. 

[11] JELENKOVIC, P., LAZAR, A. AND SEMRET, N. (1995). Multiple time scales and subexponentiality 
in MPEG video streams. Technical report CU/CTR/TR 430-95-36. Columbia University. 
http://www.ctr.columbia.edu/comet/publications. 

[12] JELENKOVIC, P. R. (1999). Subexponential loss rates in a GI/GI/1 queue with applications. 
Queueing Systems, Special Issue on Long-Tailed Distributions 33, 91-123. 

[13] KAUFMAN, J. S. (1981). Blocking in a shared resources environment. IEEE Transactions on 
Communications 29, 1474—1481. 

[14] KELLY, F. P. (1986). Blocking probabilities in large circuit-switched networks. Advances in 
Applied Probability 18, 473-505. 

[15] KELLY, F. P. (1991). Loss networks. Annals of Applied Probability 1, 319-378. 

[16] KRISHNAN, K. R. AND MEEMPAT, G. (1997). Long-range dependence in VBR video streams and 
atm traffic engineering. Performance Evaluation 30, 46-56. 

[17] LIU, L., KASHYAP, B. R. K. AND TEMPLETON, J. G. C. (1990). On the GI X /G/oo system. 
Journal of Applied Probability 27, 671-683. 

[18] LOUTH, G., MlTZENMACHER, M. AND KELLY, F. (1994). Computational complexity of loss 
networks. Theoretical Computer Science 125, 45-59. 



30 



Lu, Radovanovic 



[19] Lu, Y., RADOVANOVIC, A. AND SQUILLANTE, M. (2006). Workforce management through 
stochstic network models. Proceedings of IEEE SOLI Conference. 

[20] SEVASTYANOV, B. A. (1957). An ergodic theorem for Markov processes and its application to 
telephone systems with refusals. Theory of probability and its applications 2, 104-112. 

[21] TAKACS, L. (1980). Queues with infinitely many servers. R.A.I.R.O. Recherche Operationnelle 
14, 109-113. 

[22] WHITT, W. (1985). Blocking when service is required from several facilities simultaniously. AT&.T 
Technical Journal 64, 1807-1856. 

[23] WILKINSON, R. I. (1956). Theory of toll traffic engineering in the USA. Bell System Technical 
Journal 35, 421-513. 

[24] WOLFF, R. W. (1989). Stochastic Modeling and Theory of Queues. Prentice Hall. 

[25] ZACHARY, S . (1991). On blocking in loss networks. Advances in Applied Probability 23, 355-372. 



